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Striped high-T c superconductors such as La2-y-c C Nd y Sr cc Cu04 and La2-ccBaccCu04 near x = 1/8 
show a fascinating competition between spin and charge order and superconductivity. A theory for 
these systems therefore has to capture both the spin correlations of an antiferromagnet and the pair 
correlations of a superconductor. For this purpose we present here an effective Hartree-Fock theory 
incorporating both electron pairing with finite center-of-mass momentum and antiferromagnet ism. 
We show that this theory reproduces the key experimental features such as the formation of the 
antiferromagnetic stripe patterns at 7/8 band filling or the quasi one-dimensional electronic structure 
observed by photoemission spectroscopy. 



Unidirectional charge- and spin-density modulations 
were predicted [1] for doped transition metal oxides even 
before their experimental discovery in layered nickelates 
[2] and the rare-earth doped cuprate I^-^Sr^CuC^ [3] 
and eventually in I^-^Ba^CuC^ [4 . Stripe patterns 
emerge as a compromise between correlation driven an- 
tiferromagnetism and an optimized kinetic energy gain 
for mobile charge carriers [5]. Charge and spin stripe 
textures were indeed obtained in various approximate 
model analyses of correlated electron systems, but it has 
remained unresolved which model systems sustain stable 
groundstate solutions with stripes and superconductiv- 
ity. Here we report a pairing theory for the coexistence 
of charge and spin stripes with <i-wave superconductiv- 
ity that results from an extension of the BCS theory of 
superconductivity with an attractive pairing interaction 
for tight-binding electrons moving on a square lattice. 
Charge and spin densities and the local pairing ampli- 
tudes adjust spatially in a stripe pattern with transverse 
sign change for the antiferromagnetic (AF) order parame- 
ter. Hopping anisotropy weakens or even destroys super- 
conductivity, as observed in the low-temperature tetrag- 
onal phase of cuprate superconductors [6H2]. 

Transport experiments in the high-temperature super- 
conductor La2- x Baa;Cu04 for x = 1/8 uncovered a se- 
quence of thermal phase transitions [6[ [7J [10] . Charge- 
and spin-stripe order emerges sequentially upon cooling 
before two dimensional (2D) superconducting (SC) fluc- 
tuations set in which ultimately lead to 3D supercon- 
ductivity below 4K. These measurements provided com- 
pelling evidence for what has since been called a striped 
superconductor. The subsequently developed theory for 
the striped superconductor introduced the concept of a 
pair density wave (PDW) in which the order parameter 
for the pairing of electrons in a superconductor is spa- 
tially modulated with respect to the center of mass coor- 
dinate of the electron pair [TTJ [12]. This implies that 
Cooper pairs with finite momenta ±q form accompa- 
nied by a charge-density modulation with wavenumber 
2q [13]. The phenomenological characteristics of a pair 



density wave state with unidirectional charge modulation 
were either explored with respect to symmetry aspects 
and the nature of defects [I2j [14] or its spectral proper- 
ties [T5] . 

Although striped SC states were encountered before 
[T6UT8] . a simple microscopic model Hamiltonian which 
supports a superconducting PDW groundstate with spin 
and charge stripes has been lacking so far. Here we elab- 
orate on the existence of these solutions in an isotropic 
2D pairing Hamiltonian, characterize their real- and 
momentum-space properties, and relate them to existing 
experimental data. 

Our model is a tight-binding Hamiltonian H — Ho + 
Hi, where Ho = —^2ij^2 3 tijc\ 3 Cj 3 describes the hop- 
ping motion of free electrons on a square lattice. The 
operator cj s (ct s ) annihilates (creates) an electron on 
lattice site j with spin s =t?i; Uj are hopping ma- 
trix elements with amplitude t between nearest-neighbo, 
and t' between next-nearest neighbor sites. Here we use 
t' = — 0.4 £ for all calculations. The BCS-type attractive 
interaction 

Hl = ~^ c l c l-s c 3-sCis (1) 

is restricted to nearest-neighbor sites; V > is the pair- 
ing interaction strength. In the complete mean-field de- 
coupling scheme 

(id) (2) 
A*-A-- i 

- Vn jAi c a + 3 y 13 + V n n n oi 

we introduce the bond order parameter A^ = —V(cj^c^) 
for superconductivity and the local spin resolved densities 
rii S = {c\ s Ci S ). Using a Bogoliubov-de Gennes transfor- 
mation, the model is solved self-consistently at an elec- 
tron density 7/8 (for details on the formalism see e.g. 
Ref. |19|). The terms Vrij- s c\ s Ci S are typically not ac- 
counted for in the standard BCS theory. However, for 
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FIG. 1: Real-space characterization of the striped SC state, (a) The magnetization rrii = — exhibits AF stripes separated 
by non-magnetic anti-phase domain walls, (b) The charge density m = + m± reaches nearly one electron per site inside 
the AF stripes which are separated by single lines near quarter filling, resulting in an overall mean charge density 7/8, i.e. hole 
doping 1/8. (c) The SC bond order parameter on the horizontal (purple) and on the vertical bonds (blue). It is largest on the 
vertical bonds along the line of non-magnetic sites. Results were obtained for V = 2t on a 16x12 lattice. 



a nearest-neighbor pairing interaction they are a strong 
source for antiferromagnetism, since they provide an en- 
ergy gain —V for each AF bond, but only — V/4 for a 
bond between two non spin-polarized sites. This is the 
driving force for the formation of AF stripes in our model. 

There are two qualitatively different regimes of inter- 
action strengths: for weak V, below a critical interac- 
tion strength V c \ « 0.9 1, the only solution of the self- 
consistency equations is a homogeneous SC phase with 
d-wave symmetry and without antiferromagnetism. As 
V is increased beyond V c i, there is a sharp crossover 
into a regime where antiferromagnetism is the dominant 
order, superconductivity is suppressed and eventually 
disappears above a second critical interaction strength 
V C 2 ~ St. The characteristics of this latter regime are 
best illustrated for strong interactions V > V C 2 when Hi 
favors energetically a homogeneous AF phase for a half- 
filled band. For a filling p = 7/8 instead a configuration 
is preferred with three-legged, half-filled spin ladders and 
non-magnetic lines at an average density p = 1/2 in be- 
tween. This regularly striped solution which is unique 
for the filling p = 7/8 was indeed inferred from elastic 
neutron scattering data for the low-temperature phase of 
Lai.4 8 Ndo.4Sr .i2Cu0 4 [3]. 

If the interaction strength V is reduced below V C 2, su- 
perconductivity emerges and resides predominantly on 
the quarter-filled channels in between the AF stripes. 
Figure [T] displays the self-consistently determined mag- 
netization mi = n^ — riii ( a )> ^ ne charge density n$ = 
^it + n ii (b), and the SC order parameter (c) of the 
striped SC state for V = 2t on a 16 x 12 lattice. The 
results presented here are stable groundstate solutions ir- 
respective of the system size, provided that the selected 
geometry is commensurate with the wavelength of the 
stripes. Along the AF stripes the magnetic energy gain 
is maximized by a nearly perfect antiparallel spin align- 



ment. While stripe formation minimizes the magnetic 
energy, kinetic energy is gained by transverse fluctua- 
tions. This is the origin of the sign change in the AF 
order between neighboring stripes. 

The unidirectional character of the SC order param- 
eter is evident from its considerably smaller values on 
the bonds perpendicular to the stripe orientation (see 
Fig. [1J). The SC order parameter acquires its maximum 
value on the bonds which connect to the sites in the non- 
magnetic channels. The sign change between the SC or- 
der parameters on the horizontal and the vertical bonds 
connected to the same site verifies the d-waxe character 
of the SC order parameter. 

The groundstate solution with striped superconductiv- 
ity is not unique with respect to a sign change of A^- 
between neighboring hole-rich channels. A solution de- 
generate to the one shown in Fig. [I] exists without this 
sign change. Since in our model analysis all physical 
quantities depend on Af ■ only, these two variants of the 
striped superconductor have the same energy, provided 
Aij vanishes at the center of the AF stripes. For inter- 
action strengths close to V c i, where the AF order weak- 
ens and becomes finite also within the AF stripes, 
this degeneracy is lifted and the state without the sign 
change is favored. A similar conclusion was reached 
within a renormalized mean-field theory for a general- 
ized t-J model by Yang et a/., if the hopping amplitudes 
are anisotropic [T8] . 

A qualitative difference of the two striped SC states 
with and without sign change of the order parameter 
concerns the center-of-mass momenta q of the electron 
pairs. In the former state all pairs have either momenta 
q x = ±7r/4a or q x = ±37r/4a with q y = 0, corresponding 
to the periodicity of eight lattice sites. The state without 
sign change has in addition a finite q = 0-component for 
pairs with vanishing total momentum, i.e. (c^C-^i) ^ 0. 
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(a) momentum distribution (b) spectral density 
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FIG. 2: Momentum-space characterization of the striped su- 
perconductor (a), The momentum distribution n(k) is a su- 
perposition of incoherent states well below the Fermi energy, 
and a partially filled quasi one-dimensional band (horizontal 
bar), (b), The integrated spectral weight (see Eq. ^) cap- 
tures only the conducting states originating from the motion 
along the stripes, (c), The pair density P(k) shows the distri- 
bution of the SC pairs in momentum space, (d), The Fourier 
transformed spin density ps(k) is largest near (0, ±7r) for ver- 
tical stripes. Results were obtained for V = 2t on a 16x12 
lattice using 9x9 supercells. 



The finite center-of-mass momenta coexisting with q = 
are stabilized only by the AF stripes and vanish together 
with stripe order, thereby recovering the homogeneous 
d-wave superconductor. The state with sign change how- 
ever remains striped for sufficiently large V even in the 
absence of AF order and realizes the pure PDW state 
discussed in Ref. [T3l 

In Fig.[2]the striped superconductor is characterized in 
momentum space. The calculations were performed on a 
16 x 12 lattice with 9x9 supercells to ensure satisfactory 
momentum resolution. The momentum distribution n(k) 
in Fig. [2^i clearly exhibits the unidirectional character of 
the striped system with a horizontal bar of high occupa- 
tion probability and a diffuse region around the Brillouin 
zone (BZ) center. This diffuse background traces the 
original 2D Fermi surface of the uncorrelated electrons. 
In the absence of superconductivity the stripe order leads 
to a sharp Fermi surface with occupied states for mo- 
menta k with k y < 7r/4. The absence of discontinuites in 
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FIG. 3: Local density of states (LDOS). The colors corre- 
spond to the three distinct sites of the striped superconductor: 
sites with minimum charge density (dark blue), sites in the 
center (blue) and on the edge (light blue) of the AF stripes. 

n(k) in Fig.J^i is due to a finite energy gap in the density 
of states (see Fig. [3|. A remarkably similar momentum 
distribution has indeed been measured by Zhou et al. 
for the rare-earth doped cuprate Lai.28Ndo.6Sro.i2Cu04 
with static stripe order but no superconductivity [20] . 
For a comparison with the measured spectral weight we 
display the integrated spectral function 

^ A(k,uj)du = -- r ImG(k,k,cj)dcj (3) 

in Fig. (2}d. The lower energy cut-off at u c = —t re- 
stricts the spectral weight to the contributions of the 
non- magnetic channels (c.f. Fig. [3|. The chemical po- 
tential \i corresponds to 7/8 filling. The quasiparticle 
excitations near the Fermi level occupy the horizontal 
bar in momentum space with a strongly reduced spec- 
tral weight in the center of the BZ as observed in the 
measurements of Ref. :20j This "breach" in the spectral 
function is not captured by the physics of isolated spin 
ladders and indicates that the conducting states are not 
decoupled from the magnetic stripes. 

Similarly to Ref. [15] we define the electron-pair density 
P(k) for singlet pairing as P 2 (k) = ]T q (c_ k+q ;c kt ) 2 . 
P(k) serves as a measure at which momenta in the BZ 
electron pairs predominantly form in the superconductor. 
For the striped superconductor discussed above electron 
pairs have the finite center of mass momenta ±q where 
q y = and q x = 7r/4 or Sir/ 4 according to a stripe wave- 
length of 8 lattice constants. The pair density is ex- 
pected to be largest near the Fermi surface of the normal 
conducting system as is indeed verified in Fig. [2^. In a 
similar way we also translate the spin-stripe pattern into 
momentum space psQe) = £) q Z) s ( sc k+ qs c ks)- As snown 
in Fig. ps'(k) is strongest near (0, ±tt) for vertically 
oriented AF stripes, i.e., in those regions of the BZ where 
no electron pairs form. 

The absence of a discontinuity in n(k) is tied to the 
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FIG. 4: Doping dependence of the superconducting order 
parameter. The purple and blue curves for the maximum SC 
order parameter maxij Aij correspond to V = 2t for isotropic 



u y and anisotropic hopping t x /t y 



= 0.83, respectively. 
y with a reduced V = 



tx ty 

The grey line is the result for 

tx ty 

l.St. The dashed lines indicate an extrapolation to values of 
x — 1 where the numerical procedure does not converge to the 
discussed solution. Inset: Measured doping dependence of T c 
in La 2 - y -*Nd y Sr x Cu0 4 for y = 0.2 (Ref. 8) (red), y = 0.4 
(dark blue) , and y = (grey) (Ref. [9} • 



opening of a full gap in the DOS. The local DOS in Fig. [3] 
shows a large energy gap in the center of the AF stripes 
which is reduced on the edge of the stripe. But also on 
the non-magnetic sites, where the SC order parameter 
is strongest, there exists a small gap, because the quasi 
one-dimensionality admixes a significant extended s-wave 
component. This is in contrast to the rather 2D pure 
PDW state without antiferromagnetism, where the local 
DOS is gapless. [H CE3 US] 

The special affinity to stripe formation at the den- 
sity p = 7/8 in the cuprates is evident from the varia- 
tion of the superconducting transition temperature T c . 
The striped compounds La2- 2/ -^Nd 2y Sr :c Cu04 [8, 9] and 
La 15 / 8 Ba 1 / 8 Cu04 [3 [7] show a sharp dip in T c (x) for 
hole doping x = 1 — p = 1/8. The observed reduction 
of T c is even stronger when lattice anisotropics in the 
low-temperature tetragonal (LTT) phase grow with in- 
creasing Nd content (c.f. insert in Fig. |4|. Our model 
calculations reproduce these features as is evident from 
the doping dependence of the maximum SC order pa- 
rameter A = maxjj shown in Fig. [4] If we simulate 
the lattice anisotropy from the octahedral tilt in the LTT 
phase by introducing an anisotropy in the hopping am- 
plitudes t x ^ t y , superconductivity is weakened and the 
minimum at x = 1/8 develops into a sharp dip, which 
reaches A = for t x /t y < 0.83. For weaker pairing in- 
teraction strengths the dip at x = 1/8 develops also for 
isotropic hopping t x = t y . When x is decreased, the elec- 
tron density in the SC stripes increases, and sequentially 
at specific values of x the conducting stripes turn antifer- 
romagnetic one by one. This process is indicated by the 



dashed lines in Fig. |4j Since A is a measure of T c , the 
results for A can directly be compared to the T c data for 
La 2 - 2y - a ,Nd 2y Sr :c Cu04 012] in the inset of Fig. [I] There 
is almost no suppression of T c for the isotropic compound 
with y = 0, but a complete destruction of superconduc- 
tivity around x = 1/8 in the anisotropic Nd-doped com- 
pounds with y = 0.2 and y = 0.4. 

The presented pairing model for the coexistence of SC 
and AF stripe order reproduces the most prominent prop- 
erties of striped high-T c cuprates remarkably well. Al- 
though it was shown before that the pure PDW can be 
the groundstate of a pairing Hamiltonian in the absence 
of magnetism [13 , the issue concerning the sign change 
of the SC order parameter will only be resolved by a 
phase sensitive extension of the present calculation, e.g. 
a Josephson coupling term which has to be identified be- 
yond the Hartree-Fock decoupling scheme [XT] , Defects 
may certainly affect the stability of the stripe state con- 
siderably; the pure PDW is indeed supposed to be frag- 
ile with respect to impurities. [11] As we have verified, 
the inclusion of potential scatterers in our model shows 
that superconductivity is weakened and eventually van- 
ishes. The AF stripe order, however, is affected only 
little. Moreover, impurities can act as pinning forces for 
fluctuating stripes and thereby support the formation of 
static stripe order. 
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